Expression of circular RNAs in myelodysplastic neoplasms and their association with mutations in the splicing factor gene SF3B1.

Mutations in the splicing factor 3b subunit 1 (SF3B1) gene are frequent in myelodysplastic neoplasms (MDS). Because the splicing process is involved in the production of circular RNAs (circRNAs), we investigated the impact of SF3B1 mutations on circRNA processing. Using RNA sequencing, we measured circRNA expression in CD34+ bone marrow MDS cells. We defined circRNAs deregulated in a heterogeneous group of MDS patients and described increased circRNA formation in higher-risk MDS. We showed that the presence of SF3B1 mutations did not affect the global production of circRNAs; however, deregulation of specific circRNAs was observed. Particularly, we demonstrated that strong upregulation of circRNAs processed from the zinc finger E-box binding homeobox 1 (ZEB1) transcription factor; this upregulation was exclusive to SF3B1-mutated patients and was not observed in those with mutations in other splicing factors or other recurrently mutated genes, or with other clinical variables. Furthermore, we focused on the most upregulated ZEB1-circRNA, hsa_circ_0000228, and, by its knockdown, we demonstrated that its expression is related to mitochondrial activity. Using microRNA analyses, we proposed miR-1248 as a direct target of hsa_circ_0000228. To conclude, we demonstrated that mutated SF3B1 leads to deregulation of ZEB1-circRNAs, potentially contributing to the defects in mitochondrial metabolism observed in SF3B1-mutated MDS.

Noncoding RNAs and Their Response Predictive Value in Azacitidine-treated Patients With Myelodysplastic Syndrome and Acute Myeloid Leukemia With Myelodysplasia-related Changes.

BACKGROUND/AIM: Prediction of response to azacitidine (AZA) treatment is an important challenge in hematooncology. In addition to protein coding genes (PCGs), AZA efficiency is influenced by various noncoding RNAs (ncRNAs), including long ncRNAs (lncRNAs), circular RNAs (circRNAs), and transposable elements (TEs). MATERIALS AND METHODS: RNA sequencing was performed in patients with myelodysplastic syndromes or acute myeloid leukemia before AZA treatment to assess contribution of ncRNAs to AZA mechanisms and propose novel disease prediction biomarkers. RESULTS: Our analyses showed that lncRNAs had the strongest predictive potential. The combined set of the best predictors included 14 lncRNAs, and only four PCGs, one circRNA, and no TEs. Epigenetic regulation and recombinational repair were suggested as crucial for AZA response, and network modeling defined three deregulated lncRNAs (CTC-482H14.5, RP11-419K12.2, and RP11-736I24.4) associated with these processes. CONCLUSION: The expression of various ncRNAs can influence the effect of AZA and new ncRNA-based predictive biomarkers can be defined.

LncRNA Profiling Reveals That the Deregulation of H19, WT1-AS, TCL6, and LEF1-AS1 Is Associated with Higher-Risk Myelodysplastic Syndrome.

BACKGROUND: myelodysplastic syndrome (MDS) is a hematopoietic stem cell disorder with an incompletely known pathogenesis. Long noncoding RNAs (lncRNAs) play multiple roles in hematopoiesis and represent a new class of biomarkers and therapeutic targets, but information on their roles in MDS is limited. AIMS: here, we aimed to characterize lncRNAs deregulated in MDS that may function in disease pathogenesis. In particular, we focused on the identification of lncRNAs that could serve as novel potential biomarkers of adverse outcomes in MDS. METHODS: we performed microarray expression profiling of lncRNAs and protein-coding genes (PCGs) in the CD34+ bone marrow cells of MDS patients. Expression profiles were analyzed in relation to different aspects of the disease (i.e., diagnosis, disease subtypes, cytogenetic and mutational aberrations, and risk of progression). LncRNA-PCG networks were constructed to link deregulated lncRNAs with regulatory mechanisms associated with MDS. RESULTS: we found several lncRNAs strongly associated with disease pathogenesis (e.g., H19, WT1-AS, TCL6, LEF1-AS1, EPB41L4A-AS1, PVT1, GAS5, and ZFAS1). Of these, downregulation of LEF1-AS1 and TCL6 and upregulation of H19 and WT1-AS were associated with adverse outcomes in MDS patients. Multivariate analysis revealed that the predominant variables predictive of survival are blast count, H19 level, and TP53 mutation. Coexpression network data suggested that prognosis-related lncRNAs are predominantly related to cell adhesion and differentiation processes (H19 and WT1-AS) and mechanisms such as chromatin modification, cytokine response, and cell proliferation and death (LEF1-AS1 and TCL6). In addition, we observed that transcriptional regulation in the H19/IGF2 region is disrupted in higher-risk MDS, and discordant expression in this locus is associated with worse outcomes. CONCLUSIONS: we identified specific lncRNAs contributing to MDS pathogenesis and proposed cellular processes associated with these transcripts. Of the lncRNAs associated with patient prognosis, the level of H19 transcript might serve as a robust marker comparable to the clinical variables currently used for patient stratification.

Circulating Small Noncoding RNAs Have Specific Expression Patterns in Plasma and Extracellular Vesicles in Myelodysplastic Syndromes and Are Predictive of Patient Outcome.

Myelodysplastic syndromes (MDS) are hematopoietic stem cell disorders with large heterogeneity at the clinical and molecular levels. As diagnostic procedures shift from bone marrow biopsies towards less invasive techniques, circulating small noncoding RNAs (sncRNAs) have become of particular interest as potential novel noninvasive biomarkers of the disease. We aimed to characterize the expression profiles of circulating sncRNAs of MDS patients and to search for specific RNAs applicable as potential biomarkers. We performed small RNA-seq in paired samples of total plasma and plasma-derived extracellular vesicles (EVs) obtained from 42 patients and 17 healthy controls and analyzed the data with respect to the stage of the disease, patient survival, response to azacitidine, mutational status, and RNA editing. Significantly higher amounts of RNA material and a striking imbalance in RNA content between plasma and EVs (more than 400 significantly deregulated sncRNAs) were found in MDS patients compared to healthy controls. Moreover, the RNA content of EV cargo was more homogeneous than that of total plasma, and different RNAs were deregulated in these two types of material. Differential expression analyses identified that many hematopoiesis-related miRNAs (e.g., miR-34a, miR-125a, and miR-150) were significantly increased in MDS and that miRNAs clustered on 14q32 were specifically increased in early MDS. Only low numbers of circulating sncRNAs were significantly associated with somatic mutations in the SF3B1 or DNMT3A genes. Survival analysis defined a signature of four sncRNAs (miR-1237-3p, U33, hsa_piR_019420, and miR-548av-5p measured in EVs) as the most significantly associated with overall survival (HR = 5.866, p < 0.001). In total plasma, we identified five circulating miRNAs (miR-423-5p, miR-126-3p, miR-151a-3p, miR-125a-5p, and miR-199a-3p) whose combined expression levels could predict the response to azacitidine treatment. In conclusion, our data demonstrate that circulating sncRNAs show specific patterns in MDS and that their expression changes during disease progression, providing a rationale for the potential clinical usefulness of circulating sncRNAs in MDS prognosis. However, monitoring sncRNA levels in total plasma or in the EV fraction does not reflect one another, instead, they seem to represent distinctive snapshots of the disease and the data should be interpreted circumspectly with respect to the type of material analyzed.

Relationship between Altered miRNA Expression and DNA Methylation of the DLK1-DIO3 Region in Azacitidine-Treated Patients with Myelodysplastic Syndromes and Acute Myeloid Leukemia with Myelodysplasia-Related Changes.

The DLK1⁻DIO3 region contains a large miRNA cluster, the overexpression of which has previously been associated with myelodysplastic syndromes (MDS). To reveal whether this overexpression is epigenetically regulated, we performed an integrative analysis of miRNA/mRNA expression and DNA methylation of the regulatory sequences in the region (promoter of the MEG3 gene) in CD34+ bone marrow cells from the patients with higher-risk MDS and acute myeloid leukemia with myelodysplasia-related changes (AML-MRC), before and during hypomethylating therapy with azacytidine (AZA). Before treatment, 50% of patients showed significant miRNA/mRNA overexpression in conjunction with a diagnosis of AML-MRC. Importantly, increased level of MEG3 was associated with poor outcome. After AZA treatment, the expression levels were reduced and were closer to those seen in the healthy controls. In half of the patients, we observed significant hypermethylation in a region preceding the MEG3 gene that negatively correlated with expression. Interestingly, this hypermethylation (when found before treatment) was associated with longer progression-free survival after therapy initiation. However, neither expression nor methylation status were associated with future responsiveness to AZA treatment. In conclusion, we correlated expression and methylation changes in the DLK1⁻DIO3 region, and we propose a complex model for regulation of this region in myelodysplasia.

MicroRNA profiles as predictive markers of response to azacitidine therapy in myelodysplastic syndromes and acute myeloid leukemia.

BACKGROUND: Azacitidine (AZA) is a nucleoside analog used for treatment of myelodysplasia and the prediction of AZA responsiveness is important for the therapy management. METHODS: Using microarrays and reverse-transcription quantitative-PCR, we analyzed microRNA (miRNA) expression in bone marrow CD34+ cells of 27 patients with higher-risk myelodysplastic syndromes or acute myeloid leukemia with myelodysplasia-related changes before and during AZA treatment. RESULTS: At baseline, we found that future overall response rate was significantly higher in patients with upregulated miR-17-3p and downregulated miR-100-5p and miR-133b. Importantly, the high level of miR-100-5p at baseline was associated with shorter overall survival (HR = 4.066, P= 0.008). After AZA treatment, we observed deregulation of 30 miRNAs in responders (including downregulation of miR-10b-5p, miR-15a-5p/b-5p, miR-24-3p, and miR-148b-3p), while their levels remained unchanged in non-responders. CONCLUSIONS: Our study demonstrates that responders and non-responders have distinct miRNA patterns and that the level of specific miRNAs before therapy may predict the efficacy of AZA treatment.

Aberrant expression of the microRNA cluster in 14q32 is associated with del(5q) myelodysplastic syndrome and lenalidomide treatment.

Lenalidomide is a novel thalidomide analogue with immunomodulatory and antiangiogenic effects that has been successfully used for the treatment of low and intermediate-1 risk myelodysplastic syndromes (MDSs) with a del(5q) aberration. Because information about the influence of lenalidomide on the microRNA (miRNA) transcriptome is limited, we performed miRNA expression profiling of bone marrow CD34+ cells obtained from MDS patients with the del(5q) abnormality who had been subjected to lenalidomide treatment. To define differences in miRNA expression, we performed paired data analysis to compare the miRNA profiles of patients before and during lenalidomide treatment and those of healthy donors. The analysis showed that miRNAs clustering to the 14q32 region had a higher expression level in patient samples before treatment than in the healthy control samples, and this elevated expression was diminished following lenalidomide administration. Because some of the 14q32 miRNAs play important roles in hematopoiesis, stem cell differentiation, and apoptosis induction, the expression of this cluster may be associated with the pathophysiology of the disease.

Genome-wide miRNA profiling in myelodysplastic syndrome with del(5q) treated with lenalidomide.

OBJECTIVES: Lenalidomide is a potent drug with pleiotropic effects in patients with myelodysplastic syndrome (MDS) with deletion of the long arm of chromosome 5 [del(5q)]. We investigated its effect on regulation of microRNA (miRNA) expression profiles in del(5q) patients with MDS in vivo. METHODS: We used miRNA expression microarrays to study changes in miRNA levels in peripheral blood CD14+ monocytes collected from patients before and during lenalidomide treatment and compared them with those from healthy donors. RESULTS: Before treatment, we observed strong upregulation of pro-apoptotic miR-34a and miR-34a* that diminished during lenalidomide exposure. Upregulation of HOX-related miR-196b and erythroid-specific miR-451 seen in untreated patients remained unchanged after the treatment. At the time of hematologic response, expression of several miRNAs clustering to the 14q32 locus was reduced. Additionally, we focused more deeply on miRNAs from the 5q commonly deleted region and found that levels of miR-378 and miR-378* followed haploinsufficiency trend. CONCLUSIONS: This report describes changes in miRNA expression in del(5q) patients with MDS treated with lenalidomide, likely arising from deregulation of pathways implicated in lenalidomide action.